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Abstract. Probing the gas and dust in proto-planetary disks is central for understanding the process of planet formation. In disks 
surrounding solar type protostars, the bulk of the disk mass resides in the outer midplane, which is cold (<20 K), dense (> 10 7 
cirT 3 ) and depleted of CO. Observing the disk midplane has proved, therefore, to be a formidable challenge. Ceccarelli et al. 
(2004) detected H2D + emission in a proto-planetary disk and claimed that it probes the midplane gas. Indeed, since all heavy- 
elements bearing molecules condense out onto the grain mantles, the most abundant ions in the disk midplane are predicted to 
be Hj and its isotopomers. In this article, we carry out a theoretical study of the chemical structure of the outer midplane of 
proto-planetary disks. Using a self-consistent physical model for the flaring disk structure, we compute the abundances of 
and its deuterated forms across the disk midplane. We also provide the average column densities across the disk of , H 2 D + , 
HDJ and , and line intensities of the ground transitions of the ortho and para forms of H2D + and HDJ respectively. We discuss 
how the results depend on the cosmic ray ionization rate, dust-to-gas ratio and average grain radius, and general stellar/disk 
parameters. An important factor is the poorly understood freeze-out of N 2 molecules onto grains, which we investigate in depth. 
We finally summarize the diagnostic values of observations of the isotopomers. 
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1. Introduction 

Proto-Planetary disks are the sites of planet formation. Their 
physical, dynamical and chemical structure and evolution de- 
termine if, when, how, where and what planets form. Two very 
important parameters that are difficult to determine are the gas 
mass of the disk, and the ionization degree. Critical questions 
that need to be answered by observations are: what is the evo- 
lution of the gaseous component of the disk, in particular with 
respect to the dusty component? Is it dispersed before or af- 
ter dust coagulates into planetesimals and/or rocky planets? At 
what radius? By what process? At the same time, theory pre- 
dicts that the accretion in the disk is regulated by its ionization 
degree (Balbus & Hawley 1998, Gammie 1996). So the ques- 
tions here are: what is actually the measured ionization degree 
across the disk? What ionizes the gas? Cosmic rays and/or X- 
rays? Where do the two effects balance each other, if they do? 

Observationally answering to those questions is all but an 
easy task, especially in solar type systems. This is because the 
bulk of the disk mass resides in the outer midplane, which is 
cold and dense. As a consequence, all heavy-bearing molecules 
freeze-out onto the grain mantles, and disappear from the gas 
phase where they could be observed. So the first difficulty in 
the study of the outer disk midplane is to find probes of it. Last 
year, Ceccarelli et al. (2004) detected abundant H2D + in the 
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proto-planetary disk which surrounds the solar type protostar 
DM Tau. They claimed that H2D + probes the cold outer mid- 
plane and, in addition, its abundance is a direct measure of the 
ionization degree. The reason behind this claim is the peculiar 
chemistry of the H2D + ion. The basic idea is that in cold and 
dense gas, where CO and all heavy-bearing molecules freeze- 
out into the grain mantles, two things happen: first, only H2 and 
the ions from this molecule, namely HJ and its isotopomers, 
remain in the gas phase; second, the molecular deuteration is 
dramatically enhanced, up to having H2DVHJ larger than unity 
(Caselli et al. 2003). Therefore, Ceccarelli et al. concluded that 
a) H2D + line emission probes the gas disk midplane, and b) the 
H2D + abundance measures the ionization degree there. 

In the present article, we examine the chemistry of the 
deuterated forms of HJ in proto-planetary disks, with the goal 
of exploring, on a solid theoretical basis, the exact diagnos- 
tic value of the H2D + observations. Besides, the present study 
concerns also the other HJ isotopomers, HDJ and DJ . To ac- 
complish the goal, we develop a chemical model of the outer 
disk midplane, focused in particular on the HJ deuteration 
chemistry. The physical model that describes the disk com- 
putes self-consistently the temperature and density profiles for 
a given disk mass and star luminosity (Dullemond et al. 2001; 
Dullemond & Dominik 2004). The chemical model is based 
on what has been understood from the studies of molecular 
deuteration in pre-stellar-cores and protostars, both observa- 
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tionally (see e.g. the review by Ceccarelli 2004) and theoret- 
ically (e.g. Roberts, Millar & Herbst 2003, 2004; Walmsley, 
Flower & Pineaut des Forets 2004, Flower, Pineaut des Forets 
& Walmsley 2004). Particular emphasis is devoted to the role 
of dust grains, and the effect of dust coagulation/fragmentation 
on the disk midplane chemical structure. 

Aiming to give observable predictions, besides to provide 
the abundances of the HJ isotopomers across the disk, and the 
average column density of each species, we compute the in- 
tensities of the four ground state lines from the ortho and para 
forms of H2D + and HDJ respectively. We carry out a wide pa- 
rameter study, and explore the dependence of our results on 
three major parameters: the dust-to-gas ratio, the cosmic rays 
ionization rate, and the dust grain average sizes. Besides, we 
also discuss how the results depend on the basic properties of 
the star-disk system, namely the star luminosity, disk age, mass 
and radius. Finally, we discuss the case in which N2 disappears 
from the gas phase simultaneously with CO, which is not what 
has been observed so far, but what it would be expected based 
on the laboratory measurements of the N2 and CO binding en- 
ergies (e.g. Oberg et al. 2005). 

The article is organized as follows. In ^2] we develop our 
model. In |3]we report the model predictions of a standard 
case, and as function of the parameters of the model. In ^we 
discuss the diagnostic values of the observable quantities (line 
intensities and column densities). Finally, |5] summarizes the 
content of the article. 

2. The model 

In this section we describe the model that we developed to cal- 
culate the abundance profiles and the average column densities 
of HJ and its isotopomers across the disk, and the line emis- 
sion of the two deuterated forms of HJ which have observable 
ground state rotational transitions, H2D + and HDJ. 
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Fig. 1. The density of H2 molecules (top panel), and the local 
dust temperature in the standard disk model (bottom panel), 
described in 33.11 

cal details of the modeling procedure we refer to Dullemond et 
al (2001). A typical temperature and density profile across the 
disk (our standard case described in 33. 1> is shown in Fig^ 
The disk is flaring as can be seen from the upwards-curved 
temperature contours. We anticipate that heavy-elements bear- 
ing molecules will freeze-out, and therefore HJ deuterium frac- 
tionation will be significant at about 25 K (approximatively the 
CO condensation temperature for the involved densities) and 
below. These conditions are only fulfilled in the outer disk, ap- 
proximately outside of 20 AU. The low temperature region be- 
comes geometrically thick at large distances from the star. For 
example, at a distance of 300 AU, the 15 K contour reaches a 
height of 60 AU. Typical densities in this region are «h, = 10 7 
cm' 3 and above, up to «h 2 = 10 8 ' 5 cirT 3 in the innermost parts 
of the outer disk, near 100 AU. 



2.1. Disk structure 

For the model calculations we use a model of passively irra- 
diated hydrostatic flaring circumstellar disks (Dullemond et al 
2001; Dullemond & Dominik 2004). This model computes the 
structure (i.e. the density and temperature distribution) in a self- 
consistent way. In the center of the disk, a low-mass star is 
located with a mass of M* and a luminosity of L*. Around 
the star we distribute of material in a disk ranging from 
0.1 AU to /?disk, with a surface density powerlaw S oc r _1 , im- 
plying that the disk mass per unit radius is constant, i.e. both 
inner and outer disk contain significant amounts of mass. The 
disk contains dust at a mass fraction of /d/g which we assume to 
be fully mixed with the gas. While in reality, there probably ex- 
ists a distribution of grain sizes in the disk, we choose a single 
grains size for the present calculation. This allows us in a sim- 
ple way to study the effects of grain size. The structure of the 
disk is then computed by iterating between a 1 + 1D continuum 
radiative transfer which computes the dust temperature in the 
entire disk, and a hydrostatic equilibrium code which computes 
the vertical density and pressure distribution under the assump- 
tion that the gas and dust temperatures are equal. For the techni- 



2.2. Chemistry of deuterated HJ 

Complex chemical networks are often necessary to compute 
the abundances of molecules in interstellar environments. In 
particular in regions with a rich chemistry, it is difficult to 
predict which are the important reactions leading to and from 
a certain molecule (e.g. Semenov, Wiebe & Henning 2004). 
However, in the cold and metal-depleted regions in protostellar 
cores and circumstellar disks, the chemical network is limited 
and a much simpler treatment is possible. Also, if the chemi- 
cal timescales involved are short, a steady state solution of the 
chemistry is often appropriate. In the following we address the 
most important processes leading to the formation of deuter- 
ated HJ and develop a set of equations for the steady state so- 
lution of this network. 

In gas under standard molecular cloud conditions, HJ ions 
are formed with a rate ( by the ionization of H2 (and He) due 
to cosmic rays, and destroyed by the reactions with all neutral 
molecules and atoms in the gas, and by recombination reactions 
with electrons and grains. The reaction HJ with H2 - the most 
abundant species in molecular gas - returns HJ and therefore 
has no net effect on the HJ abundance. However, the reaction 
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with HD, the most abundant H2 isotopomer, forms H2D + : 
H3 + HD -» H 2 D + + H 2 



(1) 



with a rate coefficient k\ (see Table 1). The backward reac- 
tion is endothermic with an energy barrier of about 220 K. At 
low temperatures, it is inefficient, making the reaction^the ef- 
fective route to H2D + formation. When considering the abun- 
dances in cold molecular gas, the most important molecules 
which cause destruction of HJ are: HD, CO (the second most 
abundant molecule after H2), and N2 (see 32.31 1. Also recombi- 
nation with grains (see § !2.5t is an important process. Equating 
H+ formation and destruction rates, and ignoring the destruc- 
tion by molecules/atoms less abundant than CO and N2, the 
equilibrium abundance ratio between the HJ and H2D + is given 
by the following equation: 
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Table 1. Reactions and reaction constants important for the 
&reci*e + £co*co + £n 2 *No + k gr x„ + 2 • £-i[D] + 2&2[D]abundance of Ht and its isotopomers. The parameters y, a, 



2 • [D]*! 



(2) 



where the rate coefficients are defined in Tabled x e is the elec- 
tronic abundance, and [D] is the elemental abundance of deu- 
terium relative to H nuclei, equal to 1.5 x 1CT 5 (Lynsky 2003). 

As noticed by other authors (e.g. Caselli et al. 2003), the 
H2D + /HJ ratio, given by Equation |2 can be larger than unity 
if the gas is cold and very depleted of CO and N2. In the limit 
of very cold and completely depleted gas, it reaches fei/^ ~ 2 
(Walmsley et al. 2004). 

Just like H}, H.2D + is destroyed by reacting with molecules, 
atoms, electrons, and grains. In analogy with H}, the interac- 
tion of H_2D + with HD leads to the formation of HD} , whereas 
the reaction with H2 is endothermic and therefore suppressed 
at low temperatures. 



and p determine the rate coefficient at temperature T through 
k = y ■ (T/300K) a exp{-/3/T) and are taken from Roberts et 
al. (2004) and the UMIST database. A recent study suggests 
a lower value of the reaction rate of H} with HD (Gerlich, 
Herbst & Roueff 2002). However, the study needs confirma- 
tion, so we preferred to stick on the old value, which also of- 
fers an easier comparison with previous similar models. In the 
Table we only detail the products of those reactions produc- 
ing H} and its isotopomers. The destruction reactions involv- 
ing CO, e~ and grains may have a variety of products. In order 
to compute the rate coefficients for the reactions with CO and 
e~, we added the y parameters of all relevant reactions from the 
UMIST database. The reaction rates with grains are described 
in §0 



H 2 D + + HD 



HD} + H 2 



(3) 



As in the case of reaction Q, the reverse reaction of Q is en- 
dothermic, and therefore inhibited at low temperatures. H2D + 
is therefore destroyed by reactions with HD, CO, N2, electrons, 
and grains, and forms HD}. In the same way, the triply deuter- 
ated form of H} is formed by the successive reaction of HD} 
with HD: 



HD} + HD 



D3 



Hi 



(4) 



and destroyed by the reaction with CO, N2, grains and elec- 
trons. Now the reaction D} + HD just exchange the deuterium 
atoms, and D} is therefore the dead-end of the chain. As in 
the case of Equation [5] the abundance ratios of HDJ and D} 
with respect to H2D + and HDJ respectively can be derived by 
equating formation and destruction rates: 



2 • [D]* 2 



HD 2 + 

H 2 D + k rec2 x e + £co*co + k Nl x N2 + k gT x„ + 2 



• fe_ 2 [D] + 2k x 
(5) 



HD} 
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One can easily see that the HD}/H2D + ratio can reach unity. 
In the limiting case of very low temperatures, extreme CO and 
N2 depletion and low ionization, it reaches a value of £2/^3 ~ 
1.3. On the other hand, the abundance of D} is only limited 
by the electron abundance. In extreme conditions - cold and 
heavily CO and N2 depleted gas - D} can be the dominant 
charge carrier. 

Tabled summarizes the reactions involving H} and its iso- 
topomers, along with the rates used in this study. The reaction 
rates have been taken from Roberts et al. (2003, 2004), except 
for recombination reactions on grains which we discuss in 32.51 
In practice, the adopted chemical network is a small subset 
of the extensive chemical networks implemented by Roberts 
et al. (2003, 2004) or Walmsley et al. (2004) or Flower et al. 
(2004). Our goal here is to study the behavior of the H} and its 
isotopomers, which can indeed be described by the few equa- 
tions discussed above. This sets the limit of applicability of 
the present model to regions where CO and N2 are the only 
molecules affecting the H} chemistry. Furthermore, we did not 



£rec3*e + ^CO^CO + £n 2 *N, + £gr-% + 2£_ 3 [D] 



(6) 



treat the ortho and para forms independently, but used, when 
necessary, the results by Walmsley et al. (2004) and Flower et 
al. (2004). 
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With this reduced chemical network, we are able to repro- 
duce the results by Roberts et al. (2003, 2004), Walmsley et al. 
(2004) and Flower et al. (2004) accurately as far as the abun- 
dances of the deuterated forms of HJ are concerned. 



2.3. CO and N 2 freeze-out 

The CO and N2 abundances are obviously critical parameters 
in the deuteration of HJ. CO is known to freeze-out onto the 
grains mantles at large enough densities (> 10 5 otT 3 ) and low 
temperatures (< 25 K) (e.g. Caselli et al. 1999; Bergin et al. 
2001; Bacmann et al. 2002; Tafalla et al. 2002). Recent labo- 
ratory experiments show that the CO binding energy depends 
on the matrix where CO is embedded (e.g. Collings et al. 2003; 
Fraser et al. 2004; Oberg et al. 2005). It varies from ~ 840 K in 
CO-CO ices (average of the values measured by Collings et al. 
2003 and Oberg et al. 2005), to 885 K in a CO-N 2 ice (Oberg 
et al. 2005), and 1 180 K in a CO-H 2 ice (Collings et al. 2003; 
Fraser et al. 2004). In this study, we will adopt the value 885 
K. 

There is overwhelming observational evidence, mainly to- 
wards pre-stellar cores, that N2FF remains in the gas phase at 
larger densities than CO (e.g. Caselli et al. 1999; Bergin et al. 
2001; Tafalla et al. 2004; Pagani et al. 2004; Crapsi et al. 2005). 
Since N2H + is believed to be formed from N2, these observa- 
tions suggest that N2, which is the major nitrogen reservoir, 
freezes-out onto grains at higher densities than CO. The reason 
for that is not fully understood, because the binding energy of 
N2, measured in the laboratories, is only slightly lower that the 
CO binding energy: 790 K in a pure N2 ice and 855 K in a N2- 
CO ice (Oberg et al. 2005). However, the effect of having N2 in 
the gas phase where CO molecules disappear greatly influences 
the abundances of HJ and its isotopomers, because it causes an 
additional term of HI (H2D + and HDJ) destruction in Equation 
13 (and Eqs. |5]and|6). Therefore, we decided to adopt a semi- 
empirical approach, and to assume that the binding energy of 
N2 is 0.65 times that of CO, following the theoretical studies of 
other authors (e.g. Bergin et al. 1995, 1997). However, we will 
discuss the effect of a larger N2 binding energy, as measured in 
laboratory, in 33.41 

We treat the freeze-out of CO and N2 in a time-dependent 
way. CO molecules freeze out onto the grain mantles at a rate: 



£de P = S {na^rigj v C o 



(7) 



where we adopted a sticking coefficient S - 0.3 (Burke & 
Hollenbach 1983), and a mean grain radius of a gr . The grain 
number density n g is given by the (mass) dust-to-gas ratio /d/ g 
multiplied by the gas density n, and divided by the grain mass 
(see Eg. l llSh . A similar equation can be written for N2. 

Frozen CO and N2 molecules can be released back into the 
gas phase by thermal evaporation, and we followed the first 
order desorption kinetics approach to describe it (Hasegawa & 
Herbst 1993): 



£ ev = v Q exp[-E b /kT] 



(8) 



where Vo = 10 s is the lattice frequency, E\, is the binding 
energy per molecule of CO and N2 ice respectively, and k is the 
Boltzman constant. 



Cosmic rays also contribute to the release of CO and N2 
from the ice, with a rate of (Hasegawa & Herbst 1993): 



k CI = 9.8 x 10~ 15 



3 x lO 17 *" 1 



(9) 



The time dependent number density of gaseous CO, nco, is 
therefore given by the solution of the following equations: 



d»co 
dt 



= -kdep«CO + (^ev + ^cr) ' «C0 

«co + "co = "h 2 ■ Acq 



(10) 
(11) 



where we adopted the CO abundance Aqo observed in molec- 
ular clouds, 9.5 x 10~ 5 (Frerking et al. 1982). «™ is the num- 
ber density of frozen CO molecules. Similar equations can be 
written for N2, where the N2 abundance in molecular clouds 
is assumed to be 4 x 10 5 (Bergin et al. 1995, 1997). In our 
model, the CO and N2 abundances in the gas phase across the 
disk depend therefore on the time. They start at time=0 equal to 
their relative molecular cloud abundances. At large times, they 
are given by the equilibrium between thermal evaporation, cos- 
mic ray desorption and freeze-out onto the mantles (Leger et 
al. 1986). 

Note that we do not treat in any way the CO and N2 photo 
dissociation by the UV and/or X-rays photons emitted by the 
central star. Therefore, our model just describes the regions 
where those photons are fully shielded, i.e. the warm molec- 
ular layer and the cold midplane of sufficiently massive disks. 
As our model focuses on the HJ and its isotopomers only, it 
does not have any vocation in treating the CO and/or N2 chem- 
istry, other than computing their disappearance from the gas 
phase because of the freezing onto the dust grains. Besides, our 
model does not consider any possible grain surface chemistry 
involving CO and N2, which may transform (part of) condensed 
CO (and N2) into different, more complex molecules. 

2.4. Charge balance 

Ionization in disks can be due to cosmic rays, X rays and UV 
rays. For the current study we restrict ourselves to regions of 
the disk which are shielded from UV radiation and X rays, i.e. 
to the disk midplane. In this case, cosmic rays are the dominat- 
ing source of ionization (Semenov, Wiebe & Henning 2004). 

In regions where the heavy-elements bearing molecules 
are depleted, the positive charge is carried by the H-bearing 
species, namely HJ, its isotopomers, and H + . The negative 
charge, on the other hand, is mostly carried by free electrons. 
The fraction of electrons attached to grains remains small, in 
particular if the smallest grains have been removed by coagu- 
lation (Walmsley et al. 2004). We will therefore ignore nega- 
tively charged grains when computing the abundances of H + , 
HJ (with isotopomers), and e~. 



Similar to HJ (see §272) H + is formed by the interaction 
of cosmic rays with H2, with a rate 0.08^. Recombination of 
H + occurs mainly on grains because the direct recombination 
with electrons is too slow, and the proton exchange with HD. 
The abundance of H + is therefore independent of the electron 
density. Equating the formation and destruction rates, we have: 
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Xft+ - 



o.o8<r 

2^hd[D] + k gr x gr 



(12) 



where £ is the cosmic ray ionization rate. The recombination 
rate on grains k gr is discussed in § 12.51 

The electron abundance may be derived by equating the 
cosmic ray ionization rates with the processes removing elec- 
trons from the gas. As source term of free electrons we con- 
sider the reactions H2 + cr — > Hi + e~ and H2 + cr — > H + + e . 
As electron sink term we use the recombination of H3 and its 
isotopes with electrons, with a dissociative recombination rate 
y8, measured by McCall et al. (2003) to be 4 x 10~ 7 at 10 K. 
In addition, we assume that every recombination of H3 (and 
its isotopomers) and H + on grains remove an electron from the 
gas. This is a reasonable assumption because positively charged 
grains quickly recombine with free electrons, consistent with 
the very low abundance of positively charged grains found by 
(Walmsley et al. 2004). This leads to the following balance 
equation for the electron density: 



£(1.08>?h 2 = y6x( H ,D|+ + k gr x H + 



*(H,D)t "ff 2 'V (13) 



where X( H ,d)+ is the total abundance of and its isotopomers, 
and m is the average mass of these molecules which we set 
equal to 4.5ra H . Together with charge conservation x e - = x H + + 
xh* +*h 2 d+ +*DiH+ +*d + , a 2nd degree equation for the electron 
density results which can be easily solved. 

We can consider two limiting cases. If the electron abun- 
dance is set by the recombination of electrons with and 
its isotopes, we can expect the canonical dependence x e - °c 
-y/^/n^- If on the other hand the electron abundance is set by 
recombination of H + on grains, the electron abundance will be 
x e - oc (/n gr where n gr is the number density of grains. We will 
use this difference below to identify the main recombination 
process in the model calculations. 

2.5. Recombination on grains 

At low electron densities, collisions of ions with grains become 
an important contributer to the recombination rates. This is true 
in particular for protons, and, depending on the rate coefficients 
for dissociative recombination, also for molecular ions. At high 
densities, the charge distribution is dominated by neutral and 
singly negatively charged grains. Draine and Sutin (1987) ar- 
gue that the sticking coefficients for positive ions on both neu- 
tral and negatively charged grains should be 1. The different 
recombination rates then only depend on the Coulomb focus- 
ing factor in the collision cross section. The Coulomb focus- 
ing factor J is given by Draine & Sutin (1987) as a function 
of the grain charge, ion charge and the gas temperature. For 
a full treatment, the relative abundances of neutral and nega- 
tively charged grains would have to be calculated, as e.g. done 
by (Walmsley et al. 2004). However, in the case of really dense 
environments, in which the abundance of small grains is small 
due to the effects of dust coagulation, a simplified treatment is 
possible. Fig|2]shows the ratio J(Z gr = -l)/7(Z gl = 0). While 
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Fig. 2. Ratio J(Z gr = -l)/J(Z gI = 0) as function of the grain 
radius, for different temperatures (see text). 

for very small grain sizes and extremely low temperatures this 
ratio can be very large, in the region of parameter space stud- 
ied in this paper the factor is reasonably small. For 0.1/rni, the 
ratio is already down to 5.7. For larger grains and somewhat 
higher temperature, the ratio gets closer to one. For simplicity 
we therefore choose the following approach. We assume that a 
fraction /_ of the grains is negatively charged, and that 1 — /_ 
of the grains are neutral. The recombination rate coefficient k gr 
is then given by (Draine & Sutin 1987): 



8kT\ l/2 2 
jmiH j 



(/_/(Z = -1) + (1 - f-)J(Z = 0)) (14) 



Choosing /_ = 0.3, we expect an error in the recombination 
rates from this approach no larger than a factor of 3. 

Finally, because of the uncertainties of the grain size distri- 
bution in the midplane of a disk, we use only a single grain size 
a gr for each model. Assuming a density of the grain material of 
2.5 gr cirT 3 , the dust-to-grain ratio in number density is given 
by: 



3.2 x 10 



-12 



fd/j 
0.01/\0.1//m 



(15) 



2.6. H 2 D + and D 2 H + line intensities 

H2D + and D2H + both appear in the para and ortho forms. 
Gerlich, Herbst and Roueff (2002) have discussed in detail the 
issue, and found that at low temperatures the ortho to para ra- 
tio of H2D + is close to unity. Subsequently, Walmsley et al. 
(2004) have shown that the H2D + ortho to para ratio indeed 
varies weakly with the density in completely depleted regions, 
and it is around 0.2 at densities > 10 6 cm -3 . Since the present 
study focuses on the outer disk midplane where the densities 
are larger than 10 6 cm -3 , we do not calculate the H2D + ortho- 
to-para ratio from a chemical network, but simply assume it to 
be 0.3. Finally, Flower et al. (2004) find that the HD+ ortho- 
to-para ratio is equal to 10 for large densities and low tempera- 
tures. 
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The ortho and para forms of both H2D + and HDJ have 
ground transitions in the submillimeter to Far-IR (Tera-Hertz) 
range (Hirao & Amano 2003): 372.4 and 1370.1 GHz for the 
o-H 2 D + and p-H 2 D + respectively, and 1476.6 and 691.7 GHz 
for the o-HDJ and p-HDJ, respectively. The o-H2D + and p- 
HDJ transitions are observable with ground based telescopes 
(CSO and JCMT today, APEX and ALMA in the near future), 
whereas the other two transitions can only be observed with air- 
borne and satellite telescopes, in particular SOFIA and the up- 
coming HERSCHEL mission. In this study, we report the line 
intensity of the four transitions, computed with a non-LTE code 
which treats self-consistently the line optical depth (Ceccarelli 
et al. 2003) 1 , assuming that the disk is seen face-on. Note that, 
for the disks studied in this article, the line optical depth in the 
face-on configuration never exceeds unity though, so that the 
lines are in practice always optically thin. The collisional co- 
efficients are not well known, and we assumed that the critical 
density of all the H 2 D + and HDJ ground transitions are 10 6 
cm" 3 derived using the H 2 D + scaled collisional coefficients by 
Black et al. (1990). Since the regions where the lines origi- 
nate have relatively large densities (> 10 7 crrT3), the levels are 
mostly LTE populated, and the uncertainty on the collisional 
coefficients is, therefore, not of great importance here. The line 
fluxes are given in erg s^'cirT 2 assuming a source distance of 
140 pc. For an easy comparison with the presently available ob- 
servations, we also give the o-H 2 D + (at 372 GHz) and p-HDJ 
(at 691 GHz) line fluxes in main beam temperature, assuming 
that the observations are carried out at the CSO and JCMT tele- 
scopes respectively. 

3. Results 

In this section we discuss the results of the model calculations, 
namely the abundances of the H 2 D + , HDJ and DJ across the 
disk. One major goal of this article is to give predictions for ob- 
servations of the two deuterated forms of HJ that have ground 
rotational transitions in the submillimeter wavelength range: 
H 2 D + and D 2 H + . 

We first explore in detail our standard case, which is a disk 
with canonical values for parameters like dust-to-gas ratio, cos- 
mic ray flux and grain sizes, properties not easily derived from 
observations. In§ l3.2l we study the dependence on these param- 
eters. In § !3.3l we address system properties that are generally 
known from other observations: the total dust mass, the disk 
size, and the stellar luminosity. Finally, in § 13 .41 we discuss the 
influence of one important model assumption, the binding en- 
ergy of N 2 and CO ice. 

3.1. The standard model 

Our standard model is a disk with a dust disk mass equal to 
2 x 10~ 4 M surrounding a 0.5 L Q star whose Trff =3630 K (see 
§2.1). The dust-to-gas ratio is 1/100 in mass, and we use the 
canonical cosmic ray ionization rate £ = 3 x 10~ 17 s . We 
assume that the disk is 10 6 yr old, and extends up to 400 AU. 

1 The code has been adapted to the disk geometry. 
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Fig. 3. Abundance profiles of CO (right panel) and N 2 (left 
panel) across the disk, for our standard case. The different col- 
ors show the abundance as given by the scale bar. 



Figure^shows the density and temperature profiles across the 
disk of our standard case. 

Figure[3]shows the abundances of CO and N 2 respectively, 
across the disk. CO molecules disappear from the gas phase 
where the dust temperature is around 25 K, at a radius larger 
than about 100 AU and at heights lower than about 1/4 of the 
radius. Since the density is large in those cold regions, the drop 
in the abundance is very sharp, more than a factor 100. The 
largest CO depletion occurs between about 100 and 200 AU, 
where it is larger than 10 4 times. At larger radii it is a factor 10 
lower. N 2 , however, remains in the gas phase in a much larger 
region, and, therefore, dominates the destruction rates of HJ 
(and its isotopomers). It is only in the very outer regions, at 
a radius larger than about 300 AU that N 2 freezes-out onto the 
grains and almost totally disappears from the gas phase. The in- 
termediate region, where CO is depleted but N 2 is not, is where 
the deuteration of H J takes place in a "moderate" form, namely 
where the H 2 D + /HJ ratio is enhanced with respect to the ele- 
mental D/H abundance, but DJ is not the dominant positive 
charge carrier. 

Figure0]shows the H + , HJ isotopomers, and electron abun- 
dance profiles of the standard disk respectively. We limit the 
plot of the chemical composition to regions where the CO de- 
pletion is larger than a factor of 3. The choice of 3 is a some- 
what arbitrary. It represents when the CO abundance is similar 
to that of N 2 , or, in other words, when the role of these two 
molecules -believed to be the most abundant gaseous molecules 
at low temperatures- becomes of similar importance in the HJ 
reactions, the HJ deuteration becomes to be significant, and our 
simplified model correct. The electron abundance, is between 
3 and 30 x 10~ 10 in the outermost region, at radius > 50 AU, 
increasing going towards the upper layers of the disk. The re- 
sults are in reasonable agreement with the ionization structure 
computed by Semenov et al. (2004). In most of the disk, the 
positive charge is carried by HJ, which becomes less impor- 
tant where also the N 2 disappears from the gas phase, at radius 
>300 AU. H 2 D + is the most abundant HJ isotope, except in 
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Fig. 4. Chemical structure of the standard disk: fractional 
abundances with respect to H2 as function of the disk radius 
and height, in AU. Top left panel: H3. Top right panel: H + . 
Center left panel: H2D + . Center right panel: HDJ . Bottom left 
panel: electrons. Bottom right panel: . The color levels show 
the abundance as given by the scale bar. 

Table 2. Column densities averages across the face-on disk of 
the Hi isotopes and electrons in the standard disk model. 



Species Column Density (cm 2 ) 



H 2 D + 
HD} 



1.6 x 10 13 
5.1 x 10 12 
3.3 x 10 13 
1.0 x 10 14 



the radius >300 AU region, where D3 takes over, and become 
the most abundant positive ion. H + never plays a major role as 
charge carrier across the entire disk, even though its abundance 
is comparable to that of H2D + in outer upper levels of the disk. 

The average column densities of a face-on disk of the Hi 
isotopes and the electrons in the region where the CO depletion 
is larger than 3 are reported in Table|2] The D3 column density 
is the largest, and it is about 30% larger that of H2D + . The two 
isotopes together account for almost half of the overall positive 
charge in the disk midplane. The HDJ column density is twice 
smaller than the H2D + column density. These values are domi- 
nated by the outer region, where both CO and N2 are frozen-out 
onto the grains, and hence, becomes the dominant positive 
charge carrier. 

The line intensities of the four ortho and para H2D + and the 
HDJ ground transitions are reported in Tabled 

3.2. Ionization rate, dust-to-gas ratio and grain sizes 
as important factors for the abundances 

Important parameters of the model are the cosmic ray ioniza- 
tion rate, the dust-to-gas ratio, and the grain sizes, all of which 
are relatively ill known. One goal of this theoretical work is 
to provide predictions of observable quantities which can help 



to constrain these parameters, which all play a key role in the 
evolution of the proto-planetary disks and the eventual planet 
formation (see Introduction). 

In order to investigate the influence of the cosmic ray flux 
and of the dust-to-gas ratio, we have computed a grid of models 
where we vary these parameters within a factor 100 of the stan- 
dard value. Since the structure of the disk is fully determined by 
the dust opacity, changing the dust-to-gas ratio simply implies 
scaling the gas density while all other factors like gas tempera- 
ture and grain number density remain constant. Figure[5]shows 
the column densities as a function of the cosmic ray ionization 
rate normalized to 3 x 10~ 17 , and of the dust-to-gas ratio. It 
is immediately obvious that the electron column density varies 
as a simple powerlaw oc ^£/fd/ g °c ^£riH 2 - This corresponds 
to an electron abundance dependence x e - oc yj£/riH 2 , indicat- 
ing that cosmic ray ionization and dissociative recombination 
of electrons with the charge carrying molecules (H3 and its iso- 
topes in this case) are determining the electron abundance. This 
is consistent with the isotopes of H3 being the main charge car- 
riers. Because throughout the plot, the dust grain density n gl is 
constant (changing fy g only changes the gas density), recom- 
bination on grains would lead to a dependence x e - oc £/n Hl . 

The shape of the contours for the column densities of H2D + 
and D2H + are similar to those of e~, again indicating that these 
molecules are the dominant charge carriers. shows a similar 
behavior for very low dust-to-gas ratios (i.e. high gas densities) 
and normal or low cosmic ray ionization rates. However, for 
very low gas densities (high dust-to-gas ratios) and for very 
high ionization rates, the contours become vertical, indicating 
that the column density becomes independent of the ionization 
rate. A similar but smaller effect can also be observed for the 
column densities of H2D + and HDj at high dust-to-gas ratios. 
This is due to the fact that at low gas densities (correspond- 
ing also to low electron densities), the recombination of these 
molecules on grains starts to become important and leads to a 
faster decrease in the abundances. 

Figures[6]to0are similar to Fig.|5] but using a disk model 
with different grains sizes. We have used grain sizes of 0.01/im, 
and 1 .0/im, to study the cases of very little dust growth com- 
pared to the interstellar medium, and significant dust growth. 
Since the total dust mass in the model is constant, reducing the 
grain size corresponds to an increase in the number density of 
grains, and also in the total grain surface per cm 3 . This does fa- 





o-H 2 D + 


p-H,D+ 


o-HD+ 


p-HDj 


Transition 


1 [,o - 1 1,1 


lo.i — Oo.o 


1 1.1 _ Oo.o 


1 1,0 _ lo.i 


v (GHz) 


372.4 


1370.1 


1476.6 


691.7 


Flux erg/s/cm 2 


8.2e-18 


1.0e-16 


5.7e-17 


2.6e-18 


T m j,Av (mK km/s) 


18.8° 






4.8° 



° Note: the main beam efficiency is assumed to be 0.6 at CSO and 0.3 
at JCMT. 

Table 3. Line fluxes of the ground transitions of the ortho and 
para form of H2D + and the HDi, respectively, for the standard 
case. The velocity-integrated line intensities, expressed in main 
beam temperatures, T^Av are computed assuming observa- 
tions at CSO and JCMT of the o-H 2 D + and p-HD+ transitions 
respectively. 
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Fig. 5. Logarithm of the column density averaged over the disk 
surface of DJ (upper left panel), H2D + (lower left panel), HDJ 
(lower right panel) and electrons (upper right panel). The plots 
have been obtained for our standard case, namely with a grain 
radius of 0. 1 yum. 
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Fig. 8. Logarithm of the line fluxes as function of the dust-to- 
gas ratio and the cosmic ray ionization rate. The fluxes are in 
erg/s/cm 2 and computed for a source at 140 pc of distance. The 
four panels show the four observable transitions: o-H2D + (372 
GHz), p-H 2 D + (1370 GHz), o-HDJ (1476 GHz) and p-HD+ 
(691 GHz). 
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Fig. 6. Same as Figure|5]but for a grain radius of 0.01 fim. 
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Fig. 7. Same as Figure|5]but for a grain radius of 1 .0 yt/m. 



together, indicating a linear dependence on ionization rate and 
gas density. As discussed in § 12.41 this is an indication that H + is 
now the dominating positive charge carrier and that recombina- 
tion of H + on grains has taken over as main electron destruction 
mechanism. The reason for the significant decrease in deuter- 
ated H3 abundance compared to cases with smaller grains is 
caused by smaller CO/N2 depletions. The larger grains provide 
less grain surface for these molecules to freeze out and shift the 
abundance ratios away from the deuterated species towards HJ . 
The column densities of H 2 D + , HDJ and DJ all have a maxi- 
mum at normal cosmic ray ionization rates. If the cosmic ray 
flux is further increased, the column densities degrease again, 
because the CO/N2 depletion becomes less efficient. Also for 
smaller grains, the same effect is expected, but at much higher 
fluxes. 

Figure [8] shows the dependence of the H 2 D + and HDJ line 
fluxes on the cosmic ray ionization rate and the dust-to-gas ra- 
tio. For an easy comparison with observations obtainable at 
CSO and JCMT we also show plots with the velocity-integrated 
mean beam temperatures of the o-HDJ at 372 GHz and p-HDJ 
at 691 GHz lines respectively (Figure^- The curves of the o- 
H2D + and p-HDJ line intensities run practically parallel, so that 
the simultaneous observation of the two lines would not depend 
on the parameters of the model, but the ortho-to-para ratio of 
each species. They can therefore be "safely" used to derive the 
relative ortho-to-para ratio of H2D + with respect to HDJ . 



vor recombination on grains as a destruction route of HJ and its 
isotopes. This is indeed visible, in particular at high dust-to-gas 
ratios, the column densities of all three molecules are strongly 
decreased, and almost independent of the ionization rate. 

A different effect can be observed in the calculation for 
large grains. In this case, the column density of electrons 
greatly exceeds that of the HJ isotopes, in particular in the up- 
per left part of the diagram, at high ionization rates and high 
gas densities. The electron column density contours are closer 



3.3. The basic stellar and disk properties 

The basic stellar and disk parameters such as stellar luminosity, 
dust mass in the disk and disk size can normally derived inde- 
pendently. We therefore do not enter into a full parameter study 
of these quantities. Instead, we limit ourselves to describe the 
main trends by looking at single-parameter changes from the 
standard model. Table@]summarizes the results of these calcu- 
lations. The stellar and disk parameters mainly change the local 
densities in the disk, or the temperature. Temperature increase 
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Parameter 


o-H 2 D + 


p-H 2 D + 


o-HDJ 


p-HDJ 


L* x2 


1.21 


1.46 


1.23 


0.86 


LJ2 


0.73 


0.73 


0.73 


0.77 


M iust x 2 


1.44 


1.45 


1.46 


1.46 


M dust /2 


0.75 


0.82 


0.73 


0.60 


-Kdisk X 2 


1.64 


1.43 


1.38 


1.65 


^disk/2 


0.51 


0.73 


0.65 


0.40 



Table 4. Changes of the line fluxes emitted by the disk for 
different values of the stellar and disk parameters. The numbers 
are the fluxes of the four ground ortho and para H2D + and HD^ 
lines relative to the standard case. 



3.4. High N 2 binding energy 
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Fig. 9. Velocity-integrated line intensities, expressed as mean 
beam temperature, of the o-HDj 1 i,o - 1 1.1 (solid lines) and p- 
HD+ 1 i,o - lo,i (dashed lines) transitions for CSO and JCMT 
observations respectively. The contour plots are in mK km/s. 
Upper, center and lower panels refer to different grain radii: 
0.01, 0.1 and 1.0 yum respectively. 



influences the excitation of the lines, and may reduce depletion. 
Density changes modify the degree of ionization and again the 
depletion efficiency. The results, summarized in Table @] show 
that changes of a factor of two in any of the parameters have 
only limited effects on the line fluxes, in all cases less than a 
factor of two. The largest changes occur when the disk radius is 
changed, mostly because of the larger/smaller integration emit- 
ting region. 



As mentioned in 32. 31 there is a discrepancy between the mea- 
sured binding energy of N2 onto the grains -which is similar 
to the CO binding energy- and the observations in pre-stellar 
cores -which show that N2 molecules condense onto the grains 
at a larger density, and therefore, suggest a N2 binding energy 
lower than the CO one. In our standard model, discussed in the 
previous paragraphs, we adopted that the N2 binding energy is 
0.65 times that of CO, in agreement with the modeling of the 
pre-stellar cores observations (e.g. Bergin et al. 1995, 1997). 



j -4 -3 -2 -1 
Log(Dust-to — Gas Ratio) 

Fig. 10. Logarithm of the column density averaged over the 
disk surface of D3 (upper left panel), H2D + (lower left 
panel), HDj (lower right panel) and electrons (upper right 
panel). These plots have been obtained assuming that the N2 
binding energy is equal to that of CO. 



It is also worth noticing that observations of N2H + in a proto- 
planetary disk support indeed the lower binding energy for N2 
(Qi et al. 2003). However, for the sake of completeness, here 
we explore the influence of this assumption on the results. 

As shown in Figure [5] and previously discussed, the low 
N2 binding energy implies that N2 molecules remain in the gas 
phase in a large region where CO is frozen-out onto the grains. 
It is therefore clear that the presence or not of N2 in the gas 
phase will have a strong impact on the modeling results. 

Figure [TO] shows the column densities of H^, H2D + , HDj 
and e~ as function of the dust-to-gas ratio and the cosmic ion- 
ization ratio, similarly to Figure [5] The comparison of the two 
figures speaks for itself. If N2 has a binding energy compara- 
ble to that of CO (Figure fTOt. is the positive charge car- 
rier across most of the outer disk, in the standard case (dust-to- 
dust ratio equal to 0.01, and cosmic ray ionization rate equal to 
3 x 10~ 17 s~'), and H 2 D + and HDj are 30 times less abundant. 
On the contrary, if N2 has a reduced binding energy (Figure 
|5} the positive charge is almost equally shared between DJ, 
H2D + and H3 . This is because the N2 collisions are an impor- 
tant destroyer of the HJ , H2D + and HDj, preventing the "total 
conversion" of into . 

If the CO and N2 binding energies are similar, increasing 
and/or decreasing the dust-to-gas ratio does not change much 
the H2D + and/or HDJ column densities. In practice, increas- 
ing/decreasing the gas density with respect to the dust density 
does not change the H2D + and HDJ column densities, but just 
increases/decreases the column density -the positive charge 
carrier. This has an important observational consequence, for 
the H2D + and/or HDJ observed column densities cannot con- 
strain the dust-to-gas ratio in a large interval. However, it is 
worth noticing that, in principle, the ratio of the derived av- 
erage column densities can help to discriminate whether N2 
is present in the gas phase or not. Finally, observations of the 
N2H + in proto-planetary disks will also help to clarify whether 
N2 remains in gas phase after CO condenses (Qi et al. 2003). 
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4. Discussion 

4.1. Simplified treatment of the chemical network 

In this paper we have introduced a simplified chemical net- 
work to compute the abundances of deuterated HJ. The net- 
work treats only the four isotopes of HJ, electrons, H + , CO, 
and N2. Comparison with more extensive networks shows that 
the abundances of the {H,D}J ions can be calculated accurately 
under the conditions of significant CO and N2 depletion. The 
reason for this success is that the chemistry indeed becomes 
very simple under these conditions. CO and N2 are the last 
heavy-elements bearing molecules with significant abundances 
to leave the gas phase and to freeze out on dust grains. Together 
with dust grains, they act as the main destroyers of HJ and its 
isotopes. The simplified chemical network is, under these con- 
ditions, in fact complete in that it treats all important species 
and reactions. Under the assumption that the freeze-out of CO 
and N2 has reached its steady state, the steady state solution for 
this network can be derived directly, no time integration of the 
rate equations is necessary. The chemical network, therefore, 
lends itself for large parameter studies of cold gas. 

4.2. The main charge carrier in the disk midplane 

In the discovery paper of the H2D + line in DM Tau and TW 
Hya, Ceccarelli et al (2004) had suggested that H2D + may be 
the main charge carrier under the conditions found in a disk, 
which would allow for a direct derivation of the degree of the 
column density of positive ions by observing H2D + . With an 
H2 column density derived from dust continuum observations, 
this number directly leads to the degree of ionization, a much 
sought-after property of disk gas because of its importance for 
the magneto-rotational instability, the likely driver of disk tur- 
bulence. This conclusion has to be softened in the light of the 
present results. Our calculations show that H2D + , HDJ, DJ and 
HJ all carry similar amounts of charge, indicating that the ob- 
served column densities of H2D + and HDJ provide both an 
order-of-magnitude estimate of the degree of ionization, and 
a lower limit. It was shown earlier (Walmsley et al. 2004) that 
complete depletion in protostellar cores generally leaves H + as 
the main charge carrier. This trend, however, clearly ends at the 
densities found in disks: HJ and its isotopes contribute equally 
or dominate the charge balance. 

In this context, the binding energy of N2 is of critical im- 
portance. If N2 freezes out at the same moment as CO, then the 
disk midplane are indeed fully depleted of heavy element bear- 
ing molecules. Then, there is no important destruction route for 
DJ anymore, and the abundance of DJ sores to dominate the 
charge balance. If N2 indeed stays in the gas longer than CO, 
as suggested by observations in a proto-planetary disk (Qi et al. 
2003), H2D + and HDJ are produced in significant abundances 
and do provide a measure of the degree of ionization. 

4.3. H 2 D + and HDJ as tracers of cold gas 

One major question we aim to answer with the current study 
is: under what conditions are the observable lines of H2D + and 



HDJ tracers of the gas mass in the disk midplane? For these 
molecules to be a tracer, the observed line intensities must be 
a strong function of the dust-to-gas ratio. As we have seen in 
Figures Bland 1 101 the answer to this question depends on the 
behavior of N2. If N2 freezes-out at lower temperatures than 
CO -as supported by the observations in a proto-planetary disk 
(Qi et al. 2003)-, then the observed line intensities do depend 
strongly on the dust-to-gas ratio and can be used to measure 
it. However, under the conditions of complete depletion, i.e. if 
also N2 is completely removed from the midplane gas, the line 
intensities become insensitive to the dust-to-gas ratio. 

Fortunately, observations of N2H + in disks can discrimi- 
nate what is the reality, namely whether in the disk midplane 
CO and N2 disappear simultaneously from the gas phase or not, 
and, therefore, whether H2D + and HDJ are the main positive 
charge carriers or D J is. In the first case, observations of H2D + 
and HDJ will give a measure of the dust-to-gas ratio, once the 
cosmic ionization rate is known. This also implies that H2D + 
and HDJ measure the ionization degree in the disk midplane 
too, as assumed by Ceccarelli et al (2004). Once again, we em- 
phasize that the available observations (Qi et al. 2003) support 
this case. 

4.4. Dust coagulation 

An important result of this study is the dependence of the chem- 
ical structure of the disk midplane on the average dust grain 
size. Figure|9]shows that the line intensities of the o-H2D + and 
p-HDJ are strongly affected by the average grain sizes, if they 
are larger than 0.1 fim radius. In practice, if small dust grains 
are significantly depleted in the midplane, the line intensities 
decrease dramatically, because CO and N2 do not freeze-out 
efficiently on the reduced grain surfaces (see the discussion in 
33. 2> . and no substantial molecular deuteration takes place. The 
detection of the o-H2D + and p-HDJ at the level of the current 
instrumental sensitivity is, therefore an indication that enough 
grain surface is still available to account for the molecular de- 
pletion. Two scenarios are possible for this. First, coagulation 
might not be efficient or aggregate destruction is efficient in 
keeping the number density of small grains at reasonable val- 
ues. This is in fact in agreement with recent results of grain 
coagulation which suggest that coagulation at maximum effi- 
ciency is inconsistent with the observed opacities and shapes 
of flaring circumstellar disks (Dullemond & Dominik 2005). 
Another possibility is that the depletion of CO and N2 is very 
efficient in the beginning, before coagulation happens. Efficient 
coagulation could then form large bodies that bury the ice in- 
side and protect it from cosmic ray desorption. These consider- 
ations show that detailed models of coagulation, vertical mix- 
ing and freeze-out in combination with observations of HJ iso- 
topes can be useful diagnostics of grain coagulation. 

5. Conclusions 

We have computed the chemical structure of the midplane of 
disks surrounding solar type protostars. We have described the 
derived abundances and the resulting line intensities of the 
deuterated forms of HJ in a standard case, and for a wide range 
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of values of the dust-to-gas ratio and cosmic ray ionization rate. 
We have also discussed in detail the influence of the average 
grain sizes on the results. Finally, we have reported values also 
for different stellar luminosities, disk masses and radii, and dis- 
cussed in detail the case of the N2 freezing onto the grain man- 
tles simultaneously with CO. Our main conclusions are: 

- H3 deuteration is significant in the midplane of proto- 
planetary disks around solar type protostars. In our stan- 
dard case, the positive charge is carried by H3 in a large 
zone of the disk midplane, except at the very outer radii, 
larger than about 300 AU, where D3 takes over. H2D + is 
the most abundant H3 isotope across most of the disk, at 
radii less than about 300 AU. 

- Contrary to what was earlier assumed, H2D + is not the 
dominating positive ion in the standard case. H3 and its iso- 
topes are equally important. With increasing depletion of 
CO and N2, D3 becomes more abundant, eventually dom- 
inating as the positive charge carrier. Also, H + is only of 
minor importance. While at the pre-stellar-core densities, 
lower than 10 7 cirT 3 , H + dominates (Roberts et al. 2003, 
2004; Walmsely et al. 2004), in the disk midplane this is no 
longer the case. 

- The midplane chemical structure and the predicted line in- 
tensities of the H2D + and HDj are a strong function of the 
local cosmic ray ionization rate, which regulates the overall 
ionization degree, and the CO and N2 depletion across the 
disk. 

- The chemical structure and line intensities are also sensi- 
tive to the dust-to-gas ratio as long as depletion of heavy- 
element bearing molecules is not complete. This happens if 
the N2 binding energy is lower than the CO binding energy, 
as observed in pre-stellar-cores and possibly in proto-stellar 
disks. If, on the contrary, N2 freezes-out simultaneously 
with CO, as laboratory experiments would rather suggest, 
since the major charge carrier is D3 , increasing/decreasing 
the dust-to-gas ratio does not change appreciably the H2D + 
and HD2 column densities and/or line intensities. 

- The grain size has strong influence on the line strength ex- 
pected for H2D + and D2H + . Small grains accelerate recom- 
bination and reduce the abundances of all positive ions, and 
therefore decrease the degree of ionization. 

This article focuses entirely on theoretical predictions. A 
forthcoming paper (Dominik et al., in preparation) will analyze 
in detail the case of DM Tau, where the o-H 2 D + line at 372 
GHz has been detected, applying the model here developed to 
a practical case. Likely, with the advent of ALMA, with its 
great sensitivity and spatial resolution, observations of both the 
o-H 2 D + at 372 GHz and p-HD+ at 691 GHz will be possible 
on a routine/systematic base and with spatial resolution. Those 
observations promise to be very fruitful, and to bring unique 
information on the physical status of the midplane of the disks 
surrounding solar type protostars, likely similar to the progeni- 
tor of our own Solar System. Besides, observations of p-H 2 D + 
and o-HDj by out-of-the-atmosphere instruments will hope- 
fully measure the actual ortho-to-para ratio of these species. 
Finally, D3 would be likely the best diagnostics for the gas 



mass and ionization in low mass disks, but the only possibility 
to observe it is by absorption at ~5 /im (Ramanlal & Tennyson 
2004; Flower et al. 2004). This requires observations towards 
disks almost edge-on, but not completely, for the absorption is 
against the inner dust continuum, which must therefore be de- 
tectable. Such disks are hard to find, but new observations by, 
for example, SPITZER may well discover such sources, and 
these absorption observations may be possible in the future. 
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